arXiv:1505.08066v2 [hep-ph] 31 Oct 2015 


Scaling behaviours of the pt spectra for identified 
hadrons in pp collisions 


W. C. Zhang 

School of Physics and Information Technology, Shaanxi Normal University, Xi’an 
710119, People’s Republic of China 

E-mail: wenchao.zhang@snnu.edu.cn 

Abstract. 

We extend the scaling behaviour observed in the inclusive charged hadron transverse 
momentum (pr) distributions to the pr spectra of pions, kaons and protons produced in 
proton-proton (pp) collisions with center of mass energies (fs ) at 0.9, 2.76 and 7 TeV. 
This scaling behaviour arises when a linear transformation, p t —> p^/K, is applied on 
the pion, kaon or proton pr spectra. The scaling parameter K depends on y/s and is 
determined by a new method, the quality factor method, which does not rely on the 
shape of the scaling function. We argue that the pions, kaons and protons originate 
from different distributions of clusters which are formed by strings overlapping, and 
the scaling behaviours of these identified particles pr spectra could be understood with 
the colour string percolation model in a quantitative way simultaneously. 


PACS numbers: 13.85.Ni, 13.87.Fh 
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1. Introduction 

One of the most important observations in high energy collisions is the px spectra for 
different species of final state particles. From the spectra, we can learn a lot about the 
regularities of the particle productions. 

In many studies, searching for a scaling behaviour of the px spectra is helpful to 
reveal these regularities. In [Tj, the authors showed a scaling behaviour in the pion 
px spectra with different collision centralities at midrapidity in Au+Au collisions at 
the Relativistic Heavy Ion Collider (RHIC). This scaling behaviour was extended to 
noncentral regions in Au+Au and d+Au collisions [2], Similar scaling behaviours were 
found in the proton and anti-proton px spectra with different collision centralities at 
midrapidity in Au+Au collisions at RHIC [3J. 

Recently, we observed a scaling behaviour in the px spectra of inclusive charged 
hadrons in proton-proton (proton-antiproton) collisions at i/s = 0.9, 2.76 and 7 (0.63, 
1.8 and 1.96) TeV when a linear transformation, px —> z = px/K, was applied on these 
spectra [4J. Here the scaling parameter K depends on the collision energy y/s. This 
scaling behaviour was explained by the colour string percolation model in a qualitative 
way. In this paper, we will extend this scaling behaviour to the px spectra of pions, 
kaons and protons produced in pp collisions at 0.9, 2.76 and 7 TeV. A new method, the 
quality factor method 0E], will be adopted in the searching for the scaling parameter 
K. Compared with the method utilized in [4], this method does not rely on the shape 
of the scaling function, thus it is more robust. We will argue that the pions, kaons and 
protons originate from different distributions of clusters which are formed by strings 
overlapping, and the colour string percolation model can describe the scaling behaviours 
of the identified px spectra not only in a qualitative way but also in a quantitative way 
simultaneously. As a similar geometrical scaling behaviour was suggested when the px 
spectra of identified particles in pp collisions at 0.9, 2.76 and 7 TeV were presented in 
terms of a scaling variable r = px / Qo(px / \fs) x ^ 2 , where Q 0 = 1 GeV/c and A is around 
0.27 [7], we would like to compare the scaling behaviours presented in 2 and r. 

The organization of the paper is as follows. In section [2] we will illustrate the 
procedure to search for the scaling behaviours in the px spectra of identified particles. 
In section [3l the scaling behaviours of pions, kaons and protons will be described. Section 
[4] shows the comparison between the scaling behaviours presented in the variables 2 and 
r. I 11 section [51 we will use the colour string percolation model to explain the scaling 
behaviours of the identified particles in a quantitative way. Finally, the conclusion is 
made in section [6l 

2. Method to search for the scaling behaviours 

As done in [4], we will search for the scaling behaviours of the identified pion, kaon 
and proton p T spectra with the following steps. Let us take the pion p T spectra 
as an example, we define a scaling variable, z = px/K , and a scaled px spectra, 
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$(-) = A ■ (2vr pr) 1 d 2 N/dp T dy\ PT=Kz . Here the parameters K and A depend on the 
collision energy. By choosing proper K and A, the data points of the pion px spectra 
in pp collisions at 0.9, 2.76 and 7 TeV will migrate into one curve. As a convention, 
we usually set both K and A of the highest energy collisions to be 1. However, the px 
coverages of the pion spectra at 0.9 and 7 TeV are much smaller than the coverage at 2.76 
TeV. In order to make the scaling function <E>(z) to describe the pion spectra in the large 
Px region faithfully, we choose both K and A for the collisions at 2.76 TeV as 1. With 
different choices of K and A, the scaling functions $(*) are different. In order to get rid 
of this arbitrariness, we introduce another scaling variable, u = z/(z) = Pt/{pt), and 
the corresponding normalized scaling function T(n) = (z) 2 $({z)u)/ ff° &(z)zdz. Here 
(z) = / 0 °° z${z)zdz/ / 0 °° <&(z)zdz. The ways to search for the scaling behaviours in the 
px spectra of kaons and protons are identical to the one for pions. 

3. Scaling behaviours of identified pions, kaons and protons 

In this paper, the px spectra of identified pions, kaons and protons in pp collisions at 0.9, 
2.76 and 7 TeV were published by the ALICE collaboration [810, HU] , Here the identified 
pion, kaon and proton px spectra refer to the spectra of (7r + + 7 t“)/2, (K + + K~)/2 and 
(p+p)/ 2. The data at 2.76 TeV cover a pion px range of 0.1-20.0 GeV/c, which is much 
wider than the ones at 0.9 and 7 TeV, 0.1-2.6 GeV/c and 0.1-3.0 GeV/c. Similar results 
are obtained in the comparison among the kaon or proton px ranges at 0.9, 2.76 and 7 
TeV. As described in section [2] both the parameters K and A at 2.76 TeV are set as 1, 
thus the scaling function $(z) is nothing but the px spectra of pions, kaons or protons 
at this energy. In [4], <3>(z) is written in the Tsallis form, in which the mass effect of the 
charged hadron is ignored. In this work, since the threshold of the px range is 0.1 (0.1 
and 0.1) GeV/c at 2.76 TeV [10], which is below the mass of pions (kaons and protons) 
0.14 (0.494 and 0.938) GeV/c 2 pi], we write the scaling function <E>(z) for pions, kaons 
or protons with the modified Tsallis form 

1 

1 ~q 

( 1 ) 

where C q , q and Zo are free parameters, m is the mass of the particle species, and 1 — q 
is a measure of the non-extensivity. These free parameters are given by fitting equation 
dU) to the p T spectra of identified species at 2.76 TeV with the least y 2 s method. The 
square root of the sum of the statistical and systematic uncertainties of the data points 
has been taken into account in these fits. Table Q] lists the parameters C q , q, Zo as well 
as their uncertainties returned by the fits. The last line of this table shows the y 2 s per 
degrees of freedom (dof), named reduced y 2 s, for these fits. 


*(*) = c q 




\Jm 2 + z 2 — m 


z 0 
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Table 1 . The parameters C q , q and Zq of the scaling functions $( 2 ) for pions, kaons 
and protons. The uncertainties quoted are due to the statistical plus systematic errors 
of the data points added in quadrature. The last line shows the reduced y 2 s for the 
fits on the pt spectra of identified particles at 2.76 TeV. 



Pions 

Kaons 

Protons 

c q 

6.01T0.08 

0.214T0.002 

0.0546T0.0008 

q 

1.1416T0.0007 

1.1402T0.0004 

1.115T0.002 

z 0 (GeV/c) 

0.1308T0.0008 

0.193T0.001 

0.220T0.002 

X 2 /dof 

0.49 

0.23 

0.40 


In [4], for pp collisions, the scaling parameters K and A at 0.9 and 2.76 TeV are 
determined by fitting the scaled Tsallis distribution in equation (2) of that reference 
to the pp spectra at these two energies. The quality of the fit heavily depends on 
the parameters in the scaled Tsallis distribution, C q , q and ^o, which are fixed to the 
values obtained at 7 TeV. Thus K and A at 0.9 and 2.76 TeV rely on the shape of the 
scaling function in equation (1) of the reference. In this paper, in order to eliminate 
this reliance, we adopt the quality factor (QF) method to determine K and A for pions, 
kaons and protons at 0.9 and 7 TeV. In this method, the QF is defined in terms of a set 
of data points ( u\ v l ) 0(6] 


QF(K, A) 


E 


( v l — v l Q 2 
( u 1 — id -1 ) 2 + £ 2 


( 2 ) 


where u l = pfjIi, v l = log(A • (2'irp l T )~ 1 d 2 N l /dpfdy 1 ), the small constant £ is taken as 
0.01 and is utilized to keep the sum being finite when two points have the same u value. 
Before entering the QF formula, (u\v l ) has been rescaled so that 0 < < 1, and 

u l are ordered. Obviously, two successive data points being close in u and far in v will 
give a large contribution to the sum in equation (J2J) . As a result, a set of data points 
with a small sum (thus a large QF) is expected to lie close to a unique curve. For the 
scaling parameters at 0.9 (7) TeV, we utilize the data points at 0.9 (7) and 2.76 TeV to 
determine them. The best set of (K, A) for 0.9 (7) TeV is chosen to be the one which 
globally maximizes the QF. Table [2] tabulates K and A for pions, kaons and protons in 
pp collisions at 0.9, 2.76 and 7 TeV. In equation (Ej), the errors of data points are not 
taken into account. As described in 0, these errors can be introduced into the QF in 
the following way: 


QF(K,A) 


(p* — V 1 x f 1 p l p L 1 

( u l — u 1 - 1 ) 2 + e 2 


(3) 


where p l is the statistical uncertainty of data point i. Obviously, with this definition, 
the statistical scatter of different data points now have been taken into account in the 
QF. By maximizing the QF in equation (0), we can get another set of (Ii, A). We take 
the difference between the two sets of (A', A) returned by the maximization of the QFs 
in equations (0) and (0) as the uncertainties of ( K,A ). 
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Table 2. K and A for pions, kaons and protons at 0.9, 2.76 and 7 TeV. The errors 
quoted are taken from the difference between the two sets of (K. A) returned by the 
maximization of the QFs in equations © and ©• 



y/i (TeV) 

K 

A 


0.9 

0.93T0.04 

1.11T0.14 

Pions 

2.76 

1 

1 


7 

1.06T0.02 

0.89T0.07 


0.9 

0.913T0.007 

1.05T0.01 

Kaons 

2.76 

1 

1 


7 

1.14T0.04 

1.05T0.14 


0.9 

0.926T0.001 

1.08T0.02 

Protons 

2.76 

1 

1 


7 

1.108T0.006 

1.04T0.02 


Since the scaling parameters K and A at 0.9 and 7 TeV now are obtained with 
the QF method, we could alternatively determine the scaling function < F(^) by fitting 
equation (JH) to the combined data points at 0.9, 2.76 and 7 TeV, rather than to the 
points only at 2.76 TeV. However, the parameters C q , q and z 0 returned by this fit 
are consistent with the ones in table Q] when considering their uncertainties. Thus we 
will use these three parameters in table Q] throughout this work. With K and A in 
tabled we plot the scaled pion p T spectra at 0.9, 2.76 and 7 TeV in terms of the scaling 
variable z in the upper panel of figure [D In the log scale, we observe that all data 
points at different energies now are shifted to the same curve within error bars. This 
curve is described by the pion scaling function <l>(z) in equation (JTJ) with parameters in 
the second column of table [0 In order to see the agreement between the experimental 
data and the fitted results, we evaluate the ratios, R = (data — fitted)/data, at 0.9, 2.76 
and 7 TeV. The lower panel of figure Q] shows the R distribution as a function of z. It 
is obvious that the R values of all data points are in the range between -0.2 and 0.2, 
which implies that the data and fitted curve agreement is within 20%. This agreement 
roughly corresponds to the systematic error on R and the accuracy of the fit. If we take 
into account the systematic uncertainties on R, then the accuracy of the fit is around 
10%. Considering the fact that the data in the pion pt spectra cover about 10 orders 
of magnitude, the fit performed on the pion p T spectra is good. 

The scaling behaviours of the kaon and proton px spectra at 0.9, 2.76 and 7 TeV 
are presented in the upper panels of figures [2] and |3j As done for pions, we show the R 
distributions for kaons and protons in the lower panels of these two figures. For kaons, 
except for the last data point at 7 TeV, all the data points and the fitted curve agree 
within 20%. For protons, the data points are consistent with the curve within 20% when 
2 is smaller than 10 GeV/c. When z > 10 GeV/c, the yield and experimental errors for 
protons at 2.76 TeV have the same order. Thus the values of R are large in this high 
Pt region. 
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So far we have seen that the px spectra of pions, kaons or protons at 0.9, 2.76 and 
7 TeV indeed exhibit a scaling behaviour when these spectra are expressed in terms of 
z. The scaling function $(z) could be written as the modified Tsallis distribution in 
equation (JT]) phenomenologically. As described in [12j, the q parameter of the Tsallis 
distribution for the original px spectra of identified particles depends on the collision 
energy. In this work, the conclusion drawn from the scaling behaviours of identified 
particles px spectra is that the q parameters of the scaled (not the original) px spectra 
at 0.9 and 7 TeV are the same as the one at 2.76 TeV. Thus this conclusion is not in 
contradiction with the result from [12] . Since the scaling function <h(z) depends on the 
choice of the scaling parameters K and A at 2.76 TeV, we would like to eliminate it by 
replacing z with u — z/(z ). The (z) values for the pion, kaon and proton px spectra are 
0.442±0.007, 0.701±0.008 and 0.83±0.02 GeV/c, where the errors originate from the 
uncertainties of C q , q and zq- With the substitution of (z) and <E>(z) into T (u) defined 
in section [2] we can get the normalized scaling functions for pions, kaons and protons 
as 



where C' = ( z) 2 C q / / 0 °° &(z)zdz, q = q, u 0 = z 0 /(z) and rri = m/(z). Table [3] shows 
the values of C q , q and uq for these identihed particles. 


Pions 



z (GeV/c) 


Figure 1. Upper panel: the scaling behaviour of the pion px spectra presented in 
z for pp collisions at 0.9, 2.76 and 7 TeV. The solid curve is described by equation 
(fTj) with parameters in the second column of table [lj The data points are taken from 
HUS Ho]. Lower Panel: the R distribution as a function of z for pions. 
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Figure 2. Upper panel: the scaling behaviour of the kaon px spectra presented in 
z for pp collisions at 0.9, 2.76 and 7 TeV. The solid curve is described by equation 
(fT|) with parameters in the third column of table [l] The data points are taken from 
mmm■ Lower Panel: the R distribution as a function of z for kaons. 



Figure 3. Upper panel: the scaling behaviour of the proton px spectra presented in 
2 for pp collisions at 0.9, 2.76 and 7 TeV. The solid curve is described by equation 
dTJ) with parameters in the fourth column of table [l] The data points are taken from 
[SIHang. Lower Panel: the R distribution as a function of z for protons. The R value 
for the last data point at 7 TeV is around 0.7 and not shown in the lower panel. 


Table 3. The parameters C q , q and uq in the normalized scaling functions d/(a) for 
pions, kaons and protons. The uncertainties quoted are due to the errors of parameters 
C qi q and zq in table Q] 



Pions 

Kaons 

Protons 

c q 

3.97±0.03 

2.86±0.03 

2.28±0.03 

/ 

q 

1.1416±0.0007 

1.1402±0.0004 

1.115±0.002 

u 0 

0.296±0.009 

0.275±0.004 

0.26±0.01 


Now we would like to explore the difference among the normalized scaling functions 
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v I / Pions(w), Kaons(w) and v tp rotons (n) for pions, kaons and protons (see the upper panel 
of figure [4j) . In order to illustrate this difference clearly, we present the distribution 
of r = v h K aons(Protons)(w)/ v hpi ons (n) in the lower panel of figured In the small (large) 
p T region with u < 0.6 (u > 2), T Pions (w) is larger than T Protons (w) and ^Kaons(^), 
and the discrepancy goes down (up) with the increase of u. While in the moderate px 
region with 0.6 < u < 2, v h Pions (n) is smaller than v h Protons (u) and v h Ka ons(w). These 
differences could be validated by the comparison among the normalized moments of 
the transverse momentum distributions (pf)/(PT) n for pions, kaons and protons. This 
normalized moment is expressed in terms of T('u), (pf)/(px) n — / 0 °° u n ^(u)udu. Table[4] 
gives (pf)/{Pr) n for the identified particles with n = 2,3,4. Because of v hp ions (u) being 
above T Protons (u) and v h Ka ons(w) at small (large) u, the normalized moment of the pion 
Pt distribution is larger than the ones of the proton and kaon px distributions for small 
(large) n. Identical results are obtained from the comparison between the normalized 
moments of pions and protons produced in Au+Au collisions at y/s^N = 200 GeV [3] . 



Figure 4. Upper panel: the normalized scaling functions for pions, kaons and protons. 
Lower Panel: the distribution of the ratio between the normalized scaling functions of 
kaons (protons) and pions. 


Table 4. The normalized moments of the transverse momentum distributions 
( pf)/(PT) n for pions, kaons and protons. The uncertainties quoted are due to the 
errors of C q , q and uq in table [3j 


n 

Pions 

Kaons 

Protons 

2 

1.9T0.2 

1.7T0.1 

1.5T0.2 

3 

6.5T1.0 

4.9T0.3 

3.3T0.5 

4 

44.8T8.2 

26.8T2.2 

10.7T2.2 
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4. Comparison to the scaling behaviour presented in r 

In [7], a geometrical scaling behaviour appeared when the px spectra of identified 
particles in pp collisions at 0.9, 2.76 and 7 TeV were presented in terms of the scaling 
variable r = Pt/Qo(pt/\/s) x / 2 , where Qo = 1 GeV/c and A ~ 0.27. For the sake of 
illustrating this scaling behaviour, the author in that article chose the pion px spectra at 
7 TeV as a reference, and showed the ratio between the spectra at 7 and 0.9 (2.76) TeV 
as a function of r. In this work, in order to compare the scaling behaviours presented 
in z and r directly, we select the pion spectra at 2.76 TeV as a baseline, and plot the 
spectra at 0.9, 2.76 and 7 TeV as a function of r in the upper panel of figure El The 
pion spectra at 2.76 TeV is described by the modified Tsallis fit in equation (JT]) with z 
replaced by r. The lower panel of figure [5] shows the R distribution as a function of r. 
The data points at 0.9, 2.76 and 7 TeV agree with the fitted curve within an accuracy 
of 20%. This accuracy is comparable to the one of the scaling behaviour presented in 
z. Analogous results are obtained in the comparison between the scaling behaviours 
presented in z and r for the kaon and proton px spectra. However, the author in [7] 
showed that when r > 6, the geometrical scaling starts to break down. The data at 7 
TeV utilized in that work are preliminary and have not been published by the ALICE 
collaboration so far. In order to see whether there is a breaking in the scaling behaviour 
presented in z, more public data in the high px region at 7 TeV are needed. 



Figure 5. Upper panel: the geometrical scaling of the pion px spectra presented in 
r for pp collisions at 0.9, 2.76 and 7 TeV. The solid curve is described by the pion 
scaling function in equation m with z replaced by r. The data points are taken from 
[SI El HD] ■ Lower Panel: the R distribution as a function of r for pions. 


5. Colour String Percolation Model 

In PJH3], the author suggested the scaling behaviour presented in r could be understood 
by the gluon saturation mechanism. However, in this work, we will argue that the scaling 
behaviour presented in z for the identified particles could be explained by the colour 
string percolation (CSP) model [131 DSj in a quantitative way simultaneously. 
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In this model, colour strings are stretched between the partons of protons in pp 
collisions. These strings then decay into new ones with the emission of qq pairs. 
Identified hadrons are produced through the hadronization of these new strings. The 
strings are viewed as small areas of Si = nrf with r 0 ~ 0.2 fm in the transverse plane. 
The number of strings grows with the increase of the collision energy. When there are 
n strings, they start to overlap and form clusters with transverse areas of S n . The 
mean transverse momentum squared (pf)ni of identified particles produced by a cluster 
can be written as (pf)ni = \j'nSi/S n (px)u, where {pf)u is the mean pf of identified 
particles produced by a single string, and nSi/S n is the degree of string overlap. If 
strings just get in touch with each other, then S n = nSi and nS\/S n = 1. If strings 
maximally overlap with each other, then S n = S i and nSi/S n = n with n > 1. The 
Pt spectra of identified particles produced in pp collisions is deemed as a superposition 
of the px spectra produced by each cluster, f(x,px), weighted with the distribution of 
the cluster’s size, W(x). Here the cluster size x is equivalent to the inverse of ( pf) n i. 
In [IT] and |KJ, the fragmentation function for the cluster is chosen as the Schwinger 
formula f(x,px) = exp (—pfx) [16]. However, this function only describes the spectra 
well at the soft px region. In order to describe the spectra at the high px region as well, 
the Schwinger formula should be replaced by f(x,px) = exp(—y / 2xpr) [HUH]- W{x) 
is supposed to be a gamma distribution, W{x) = p^y( 7 x) fe ^ 1 exp(—qx), where 7 and k 
are free parameters. Thus the px distribution of identified particles in pp collisions is 



(5) 


where C is a normalization constant for the total number of clusters formed for identified 
particles before hadronization. In order to see whether the CSP model can describe the 
scaling behaviour presented in z, we fit equation (j5]) to the px spectra of identified 
particles at 2.76 TeV. The parameters C, r and k returned by the fits as well as the 


reduced y 2 s are tabulated in table [5j As shown in the upper panel of figure [6j the data 


points for pions agree with the CSP fit well in log scale. The accuracy of this agreement 
is around 20%, which can be seen from the R distribution in the lower panel of the 


figure. Similar results are obtained for kaons and protons. As the k parameters are 


different for the identified particles, pions, kaons and protons originate from different 
distributions of clusters. 
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Table 5. The parameters C, r and k of the CSP fits on the spectra of pions, kaons and 
protons at 2.76 TeV. The uncertainties quoted are due to the statistical plus systematic 
errors of the data points added in quadrature. The last line shows the reduced y 2 s for 
the fits. 



Pions 

Kaons 

Protons 

c 

9.1T0.2 

0.53T0.01 

0.160T0.007 

7 

0.140T0.003 

0.342T0.008 

0.91T0.07 

k 

3.15T0.02 

3.33T0.02 

5.0T0.2 

x 2 / dof 

1.53 

0.91 

1.91 



Figure 6. Upper panel: the scaling behaviour of the pion pt spectra presented in z 
for pp collisions at 0.9, 2.76 and 7 TeV. The solid curve is described by the CSP fit 
in equation (0 with parameters in the second column of table [5] The data points are 
taken from la Ga mg. Lower Panel: the R distribution as a function of z for pions. 

Now we would like to seek for the reason why the CSP model can describe the 
scaling behaviour of the identified pp spectra. Under the transformation x —> x' = Ax, 
7 — > y — 7 /A and pp —>■ pf = Pt/V A, both W(x) and f(x) are invariant. So the 
Pt spectra of identified particles in equation (J5]) are also invariant. This invariance is 
the scaling behaviour we are looking for. Comparing the p' T transformation in the CSP 
model pf — > p' T \J A with the one used to search for the scaling behaviour p T —y pr/K, 
we know that K is proportional to l/VX As described in |jT5j, HE], A = (. S n /nS \) 1 ^ 2 , 
where the average is taken over all clusters for identified particles, thus K should be 
proportional to (nS'i/S'n) 1 / 4 . Since the degree of string overlap nS\/S n grows with 
energy, the scaling parameter K should also increase with energy. That’s indeed what 
we see for the identified particles in the third column of table [ 2 ] Therefore the CSP 
model can explain the scaling behaviours presented in z for pions, kaons and protons 
separately. However, the rate of the K values increasing with energy for pions is different 
with the one for kaons or protons. This could also be explained by the CSP model. As 
K = ( p T )/ (z ), and the values of (z) are the same for pions (kaons or protons) at 0.9, 
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2.76 and 7 TeV, the ratio between the values of K at different energies should be equal 
to the ratio between the values of (pfi) at different energies, (pfi) is evaluated in terms 
of the CSP model as |7j 

, v /o°° f™W(x)f(x,p T )pfdxdp T 

[PT) /o°° /o°° W(x)f(x, p T )pTdxdp T • 1 J 

The cluster fragmentation functions f(x,pr) are the same for pions (kaons or protons) 
at different energies. However, the cluster’s size distributions W{x) are different for 
different particles, which can be seen from the comparison among the k values for pions, 
kaons and protons in table [5j As a result, for different species of particles, the values 
of K change in a different rate with energy. In a summary, the CSP model could 
be utilized to explain the scaling behaviours presented in z for the identified particles 
simultaneously. 

6. Conclusions 

In this paper, we have extended the scaling behaviour in the inclusive charged hadron 
Pt spectra to the spectra of identified particles at 0.9, 2.76 and 7 TeV. This scaling 
behaviour is exhibited when pr is replaced by z — Pt/K. The scaling parameter K is 
determined by the quality factor method, which does not rely on the shape of the scaling 
function. We have argued that the pions, kaons and protons originate from different 
distributions of clusters formed by strings overlapping, and the scaling behaviours of 
these identified particles could be explained by the colour string percolation model in a 
quantitative way at the same time. 
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